Disordered one-dimensional Bose-Fermi mixtures: The Bose-Fermi glass 
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We analyze an interacting Bose-Fermi mixture in a ID disordered potential using a combination 
of renormalization group and variational methods. We obtain the complete phase diagram in the 
incommensurate case as a function of bosonic and inter-species interaction strengths, in the weak 
disorder limit. We find that the system is characterized by several phase transitions between su- 
perfluid and various glassy insulating states, including a new Bose-Fermi glass phase, where both 
species are coupled and localized. We show that the dynamical structure factor, as measured through 
Bragg scattering experiments, can distinguish between the various localized phases and probe their 
dynamics. 

PACS numbers: 67.60.Fp, 67.85.Pq, 71.10.Pm, 71.23.-k, 71.55.Jv 



Introduction. Spectacular developments of experimen- 
tal techniques in the field of ultra-cold atoms have opened 
a new way in the study of strongly correlated systems 
[l||. Experimentalists are now able to create optical lat- 
tices with highly tunable parameters, and realize vari- 
ous models and quantum phase transitions. In Ref. [3, 
e.g., superfluid bosons were driven to a Mott insulating 
state, by tuning the on-site interactions of the Hubbard 
model, and the incompressibility of this state has also 
been demonstrated [3J. The Mott transition was also re- 
alized in Fermi systems '4], where other interesting phe- 
nomena such as pairing with spin imbalance [5, J3] or the 
BCS-BEC crossover [1] were also investigated. 

Having more and more control over conventional sys- 
tems, experimentalists and theorists now turn to the 
study of more complicated ones. On the one hand, a 
lot of attention is devoted to multicomponent systems 
such as Bose-Fermi mixtures [THllj. or three component 
systems 12J, which offer the possibility of realizing new 
phases of matter, such as supersolids, color superconduc- 
tivity |13|, or 'baryonic' phases [1J|. On the other hand, 
the creation of disorder using speckle lasers or incommen- 
surate laser beams in trapped systems paved the road to 
create quantum glasses and Andersoir insulators in cold 
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atomic systems 15|, |l6| 



It is well known that impurities can lead to the local- 
ization of electronic wave functions. While interactions 
are known to be importairt in the localized phase, taking 
them into account is quite a difficult task. In one dimen- 
sion, this question has been addressed both for fermions 
|17| and for bosons [17, 18]. For spinless fermions, suffi- 
ciently attractive interactions trigger a transition from 
the Anderson insulating phase to a p-wave supercon- 
ducting phase. For bosons, a disorder-induced insulating 
phase appears for sufficiently strong interactions: nick- 
named Bose glass, it has no gap, it is insulating, and 
is compressible [18|. The Bose glass disappears at suf- 
ficiently small repulsive interactions, where the system 
enters a superfluid phase. On the other hand, for eveir 



FIG. 1: Phase diagram of a ID *''Rb-'*''K Bose-Fermi mix- 
ture, in the weak disorder limit. We consider an array of 
tubes, created with lasers of wavelength A = 1064 nm, which 
corresponds to a (2D) lattice constant d — A/2. We take a;,;, = 
lOOao as Bose-Bose scattering length, while the Bose-Fermi 
scattering length, Obf, is tuned using a Feshbach resonance. 
The one. dimensional densities are chosen to be pfd — 0.3 
and pbd — 0.1. The recoil energy is -E_r,(, — h^\~^ /{2miib), 
while V±b is the transverse confining potential creating the 
ID tubes. The Bose-Bose interaction Ub increases with V±b 
[301 . Four phases and the region of instability of the Lut- 
tinger liquid theory are shown. BFG: Bose-Fermi glass, BFG* 
(BFG with an extremely large bosonic localization length), 
AG+SFB: Anderson Glass + Superfiuid Bosons, LL: Lut- 
tinger liquid (above the instability region) . 



smaller interactions bosons should undergo Anderson lo- 
calization and condense in the lowest eigenstate [ITI . Il8| . 
The Bose glass [19] and the Anderson glass [15| phases 
were now observed experimentally. 

In the present paper we wish study, how the most basic 



two-component system, the Bose-Fermi mixture (BFM) 
is affected by tfie presence of disorder. Tliis is not only 
important in view of the current experimental interest 
|7h11|. but it is a challenging fundamental question on 
its own, too, especially in the limit of strong interac- 
tions. In particular, one may wonder whether disorder 
"decouples" both species or whether the localization of 
one species triggers the localization of the other species 
because of interactions. This issue has been discussed 
in a general context in Ref. [20|, where a variety of ex- 
otic situations were considered. In this letter we focus 
on a disordered one-dimensional (ID) BFM, where inter- 
actions are particularly relevant. One-dimensional BFMs 
in the absence of disorder have been investigated by a va- 
riety of analytical 2lL l22l| and numerical [23| techniques. 
Here, we follow Refs. [ITI . ll8J and use bosonization to 
extract the low-energy properties of the BFM in a disor- 
dered potential. We focus on the case where the bosonic 
and fermionic densities are not commensurate either with 
the lattice constant or with each other. Our main results 
are summarized for a ^^Rb-'^'^K mixture in Fig. [TJ We 
find two different glassy phases as well as a Luttinger liq- 
uid phase and a region where the Luttinger liquid theory 
breaks down (see caption). 

Model. We consider a ID mixture of bosons and spin- 
less fermions described by the following Hamiltonian, 
H = Hb + Hf+Hbf+ iJdis, where 



Hn 



dx 



-^iix)^M^) + ^:plix) 



Hbf = Ubf / dx Spb{x)Spf{x), 



Ha;s = 



dx 



(Vb{x)pb{x) + Vf{x)pf{x) 



(1) 
(2) 



Here ipa (a = b, /) denote the Bose/Fermi fields, pa{x) = 
iljl^{x)ilja{x) their density, and nia the mass of the 
bosons/fermions. The symbol : • • • : stands for normal 
ordering. The interaction between bosons must be re- 
pulsive for stability reasons, Ub > 0, while for spinless 
fermions with a local interaction, Uf = 0. Bosons and 
fermions can have either repulsive {Ubf > 0) or attractive 
{Ubf < 0) interactions. The random potentials, Vf/b{x) 
describe the effects of disorder. For simplicity, here we 
assume them to be Gaussian with zero mean values, and 



spatially uncorrelated, i.e., Va{x)Va{x') = Da5{x — x') 
[3l|. Experimentally, there is only one disorder poten- 
tial, for instance a laser speckle. As a result, Vf and Vb 
are correlated |32l |. Nevertheless, pb being incommensu- 
rate with pf, we find that for long wavelength excitations 
Vb and Vf act as if they were independent [28j. In the 
following, we therefore assume Vf{x)Vb{x') = 0. 

In ID, long wavelength excitations are essentially den- 
sity fluctuations, and the densities Qf ib and the fields 
'0//b can be "bosonized", i.e., represented in terms of 



bosonic phases, 'i'flh 
pie quadratic forms 



[24]. Then Hf and Hb take on sim- 



Ha — 



27r 



dx 



Ka {dxOa{x)f + -TT- {dx(l)a{x)f 
Ka 



, (3) 



with (j)a{x), 9a{x) denoting density and phase quantum 
fields, obeying [(j)a{x'),dxOa{x)\ — iT:S{x ~ x'). The pref- 
actor Vb can be viewed as the sound velocity for Bo- 
goliubov phonons in a quasi-condensate, while Vf is the 
Fermi velocity. The dimensionless Luttinger parameters, 
Kf and Kb, characterize the strength of interactions [25j: 
In our case of non-interacting fermions, Kf = 1. For 
bosons, Kb ^ 1 in the hardcore limit, Kb > 1 for softer 
repulsive interactions, and Kb -^ oo for non-interacting 
bosons. Keeping in mind that gb is incommensurate with 
Qf, backward scattering of bosons on fermions is irrele- 
vant, and at low energies Hbf can be approximated as 



Hbf - 



^2 



dx dx4>b{x)dx(t>f{x) 



(4) 



In the same way, disorder induces backward scattering 
and forward scattering. However, the forward scattering 
can be gauged away (without affecting current-current 
and superfluid correlations), and only backward scatter- 
ing responsible for localization remains [17, 28iJ, 



TTOt 



da;L„^„(x)e-*2'^°(^)-t-/i.c. 



(5) 



with S,c{x)Ca{x') = Dc 5{x - x'). 

Renormalization group (RG). To perform the RG anal- 
ysis, we use the replica trick. We introduce n replicas of 
the system, average over disorder, and then take n — > 0. 
The replicated action then reads S'rcp = 5*0 + S''^^, with 



'^-T^X^^Jd-d- 



a=l a=f,b 



-l(5.C)'+^a(5.C)' 



J2^ f dxdT d^rfdxcl^t, 



(6) 



-^Z = -DfP} J2 I dxdrdr' cos{2<j)}{x, r) - 20^(x, r')) - Dbpl ^ f dxdrdr' cos(20;j(x, r) - 2(f>l{x, r')). (7) 

a,b a.b 



r 



Integrating out high momentum degrees of freedom we 



find the following flow equations to lowest order, 
dDa 



dl 



= {3-X^)D^{l), 



-AK^v^) = -2K^{l)va.{l)D^{l) C^ 
dl 



(8) 
(9) 



while -^{Ka/va) = 0. Here we have defined dimen- 

2 K^ 
sionless variables for the disorder, Da — — : T'^a, 

ttAq vi 

and parametrized the high-energy cut-ofF as A = Aoe~'. 
The Ca are numbers coming from our RG scheme. The 
anomalous dimensions, Xjii,, of the disorder operators 
can be obtained from the diagonalization of Sq and read: 



Xf 






(10) 



where t — Vf/vf, and g — Uhf /t^^/ KfKi,/ {vfV})) is a di- 
mensionless parameter. The dimension Xi, is obtained 
by changing t — >■ 1/i and Kf into Kb in Eq (fTO|) . 

For uncoupled species {g = 0), Xf — 2Kf and Xb — 
2Kb. Thus spinless fermions (bosons) are localized when 
Kf < 3/2 {Kb < 3/2) [l3|. At weak disorder we can 
neglect the feed-back, Eq. ©, and follow, e.g. Ref. |25, 
to relate Kf, Kb, t, g to experimental parameters, and 
construct the phase diagram of Fig. [TJ The thick solid 
blue line denotes Xf = 3, while the dashed red line in- 
dicates Xb — 'i- Naively, one would expect these lines 
to separate three phases: In the Bose-Fermi glass phase 
(BFG) Xb < 3 and Xf < 3, both species are pinned by 
disorder, and fermionic and bosonic excitations are local- 
ized over the localization lengths, ^/ and ^b, respectively. 
In the Luttinger liquid (LL) phase Xb > 3 and Xf > 3, 
disorder is thus irrelevant, and both species are super- 
fluid. In the regime, Xb > 3 and Xf < 3 one would 
naively predict a phase with localized fermions forming 
an Anderson glass (AG) and superfluid bosons (SEE). In 
this regime, however, some care needs be taken: since 
fermions are localized, fermionic density fluctuations be- 
come "gapped" at a length scale A ~ 1/C/- Below this 
scale bosons interact with their bare interaction. Thus 
the RG equations become 



dlogPb 
dl 



i-Xb 

3 -2Kb 



if A > 1/^f 
if A<l/e/. 



(11) 



As a consequence, the region Xb > 3 and Xf < 3 is 
divided into two phases. In the region, Xf < 3 and 
Kb > 3/2, we obtain an AG+SFB phase. However, in 
the region Xf < 3 and Kb < 3/2 bosons are ultimately 
localized on a very large length scale, and we find a Bose- 
Fermi glass phase (BEG*). 

To support the RG picture and to describe the 
"gapped" phases, we made use of the Gaussian varia- 
tional method (GVM) in replica space, as introduced in 
[26'| and used in Ref. [27j to treat interacting ID disor- 
dered media. This method is able to capture the localized 
compressible phases. To start with, we rewrite the action 
So in Fourier space as 



5" 



1 1 
2jL 



E 



(l)l{q,iuJn){GQ')%{q,iuJnWp{-q,~iUr, 



where a, 13= f, b while Latin indices run from 1 to n, the 
number of replicas. The free propagator is a 2n x 2n ma- 
trix given by {GQ'^)'^''a{q,iUJn)=6ab/{T^Ka){ujl/Va+Vaq'^) 

and {G^^)j'i{q,iuJn) = iGa')t^fiq,iLO„) = SabUbfq^ h^ 
The idea of the GVM is to replace the complicated action 
S in replica space by its best Gaussian approximation, 

and CT° o the self- 



5g, with {G-^)fa = {G^^)% - a 



energy. The optimal G can then be obtained by minimiz- 
ing the variational free energy, Fya^=Fa + {S — Sg)g/(3 
with respect to G'^o. 
Similar to Ref. 



27| . we find that the phase with lo- 



calized fermions and superfluid bosons is well described 
by assuming one step replica symmetry breaking (IRSB) 
in the fermionic sector [28]. We obtain a fermion mass 



-/ 



c 



/ 



S/ 



D, 



yT^5^A 



A^ 



(12) 



in agreement with the RG, and the solution ceases to 
exist when Xf > 3. An RG calculation with respect to 
the latter variational action {Sg with ^f finite) provides 
a modified RG equation which indeed smoothly interpo- 
lates between the limiting cases in Eq. ([Tl]) . and can be 
used to compute the localization length, ^b- In the in- 
termediate region BFG* we find an extremely large, but 
finite ^h (see Fig. [2]). Similarly, we can compute the su- 
perfluid stiffness V = lim liTauj'^{(pb{q,^)4'b{^q, — w)) in 

the AG+SFB phase using our variational solution. We 
flnd that V > only if Kb > 3/2, indicating a phase 
transition towards a Bose Glass phase. This BFG* phase 
for Kb < 3/2 appears in the variational problem as a a 
level 2 replica-symmetry breaking solution (RSB), with 
the replica symmetry also broken in the bosonic sector. 
However, ^b being extremely large in the BFG* region, it 
may appear as a superfluid phase in a flnite system [28J. 
The nature of the transition between the BFG and the 
LL is quite subtle and will be detailed in J28| . 

Observables. The dynamical response of the 

Bose-Fermi mixture varies from one phase to an- 
other. Such a response can be probed through Bragg 
scattering measurements (see e.g. [29|), giving ac- 
cess to the dynamical structure factor Sb/f{q,i^) = 
/ dtdx e*«^-'"* {Pb/f{x, t)pb/f{0, 0)). The variational ap- 
proach enables us to compute Sb/f (?, w) , at least for small 
values of q. Fig. [2] shows the inelastic part of Sf{q,uj) 
computed at two points of the phase diagram. The 
elastic part that comes from the gauge transformation 
used to absorb forward scattering would appear as a 
Dirac delta at lu=0. The main plot shows the fermion 
structure factor in the BFG phase. It exhibits a two 
peak structure with peak locations cj w Wf,(g^ -I- T,b)^^'^ 
and uj w Vf{q'^ + E/)^/^. The first peak is a bosonic 
feature, indicating that the Bose-Fermi glass is indeed 
a coupled localized phase. The main effect of disor- 
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FIG. 2: Main plot: femiion structure factor Sf{q — 0.2Ao,tj) 
in the BFG phase as computed by GVM. Bose-Fermi scatter- 
ing length is atf — — 35ao. The UV cut-off is Aq = pb- We 
have taken V±^b = 4^0Er,i,. Other parameters are the ones 
of Fig. [T] Bottom inset is a zoom around the first peak at 
oj = Vb{q'^ + 'SbY , a signature of the coupled nature of the 
Bose-Fermi glass. The additional dashed curve is the struc- 
ture factor in the AG+SFB phase (VI, 6 = 13Eji,b)- Top inset 
shows Efc = ^(^ , as obtained from the RG fiow (circles) and 
the 2-step RSB variational calculation (triangles). Vj_,b (i.e. 
Kb) is fixed and the boundaries obtained from the first RG 
analysis are represented. One can see that in region BFG*, 
the mass is non zero but extremely small. Below abf — — 95ao 
{Xf =3), Efc is identically zero. 



der is to introduce new energy scales in the structure 



factor, VfY^ 



1/2 
/ 



and ffeSj . 



Its other effect is to intro- 



duce a linear frequency dependence for small lv |27l . 
Notice that in the localized phases there is no hard 
gap in the excitation spectrum. The bottom inset is a 
zoom around the bosonic peak. Its counterpart in the 
AG+SFB phase is plotted there (for clarity it has been 
left out in the main plot). Note that in the AG+SFB Efc 
is zero and the peak at a; ~ v^q is much sharper, as it 
would be in a simple Luttinger liquid. Finally, in the LL 
phase, where both species are free, one would simply find 
two very sharp peaks, corresponding to the two sound 
modes of the mixture, at a; = V-i^q and uj = V-q, with 

1/2 



V 



v\ 



2^,2 



The transition from the AG 



fv'jvl 



21|,|22 



-SFB state to the BFG 
phase can easily be detected through time of flight 
(TOF) measurements. In the AG+SFB phase one has 



{i,,{xWb{x')) 



l/\x ~ a;'|^/^^% with Kl the renormal- 
ized Luttinger parameter of the bosons. As a result, a 
bosonic coherence peak is predicted with a power law de- 
pendence, nb{R) ~ l/R^~^/^^>' , with nt, the density of 
bosons as measured in the TOF experiment, and R the 



expansion along the ID tubes. 

Conclusion To summarize, we have established the 
phase diagram of an interacting Bose-Fermi mixture in 
the presence of uncorrelated disorder. For relevant exper- 
imental parameters, in the case of a ^^Rb-'**'K mixture, 
we have found three different phases, including a Bose- 
Fermi glass where both species are coupled and localized. 
These new phases can be detected by Bragg scattering 
and time of flight measurements. 
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